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The Matrix Interpretive System devised by Wilson and 
modified by Cantin has been transformed into an interactive 
program operable under a time sharing system. Several new 
commands have been added to extend the usefulness of the code 
and give it the capability to offline moods write and punch 
data, to plot graphs and contour maps and to integrate simple 


polynomial expressions. 
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i. ~ INTRODUCTION 


A. GENERAL DESCRIPTION 

Many engineering problems in such fields as steady-state 
temperature distribution [Ref. 1], stress in elastic structures 
[Ref. 2], and fluid flows [Ref. 3] can be formulated as problems 
of numerical linear algebra. Until the advent of modern 
digital computers, all but the simplest of these problems had 
to be reserved to specialists for solutions. Even problems of 
moderate size required a considerable computational effort. 
The computer changed all of this. All of the simple matrix 
operations like multiplication, transposition, solution of eigen 
systems can be coded to solve individual problems. The computer 
program described in this report integrates a number of stand- 
ard subprograms of linear algebra [Refs. 4,5,6,7] into a 
problem oriented computer language. This program, heretofore 
called: Matrix Interpretive System was coded in FORTRAN IV. 
The resulting matrix operations sub-language has a syntax and 
rules of its own as well as a certain number of diagnostic 
messages that can be automatically returned to the user of the 
language. The basic specification for the design of the 
language structure was that it had to be simple enough to be 
used by persons unfamiliar with computer programming. As 
written, the language could be used to solve stacks of problems 
sequentially; if an error is found in one of the problems a 
message is printed and a solution for the next problem is 
attempted. To be effective in a conversational mode at a 


remote station, the program had to be heavily modified so as 
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to leave effective control of the flow of operations to the 
user. This version of the code is called: Interactive Matrix 
Interpretive System. 

The Interactive Matrix Interpetive System is a general 
program written in FORTRAN IV [Ref. 81, it is a problem oriented 
Sub-language. The entire program is coded in double precision 
arithmetic. The only subroutines which are used but are not 
contained in the program are a square root and a natural log 


subroutine. The program uses 160,000 bytes of core storage. 


B. HISTORICAL BACKGROUND 

The original program was developed by Professor E, L. Wilson 
©: the University of California in April 1963 for an IBM 7094. 
It was later adapted by Professor G. Cantin of the Naval 
Postgraduate School for a CDC 1604 in February 1966 and for 


the IBM 360/67 in August 1968. 


Cc. OBJECTIVES 
The objectives of the author's work have been to: 
ie transform the program into an interactive system 


operable under time sharing from a remote typewriter station; 


is add operations to extend the capabilities of the 
program; 
or create a "USER'S MANUAL" to facilitate the use of 


the program. 
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II. GENERAL PROGRAM OPERATION 

The Interactive Matrix Interpretive System will be referred 
—O aS IMIS in this and following sections of this report. IMIS 
is written in FORTRAN IV language for the IBM 360/67 time 
sharing system available at the Naval Postgraduate School 
(CP/CMS). A different version called DSMIS is used for batch 
processing. All operations are performed in double precision 
arithmetic. Matrices are stored internally in a single 
dimensional array, A, that is 10,000 double words in size. 
Matrix names are stored in an array called NAME. Arrays, NROW 
and NCOL, are used to store the number of rows and columns of 
each matrix. 

The following is a brief description of the function and 
operation of most of the subroutines of IMIS. Some of the 
subroutines are straightforward and will not be discussed. A 


flow diagram of the program is presented in Figure l. 


A. MAIN PROGRAM AND SERVICE SUBROUTINES 
i. Main Program 
The Main program is used only to initialize internal 


variables and to call the main service subroutine, . SMIS1. 


2s Main Working Subroutine (SMIS1) 


SMIS1 is the general working subroutine. Some of the 
operations such as duplication of a matrix and envelope of 
a matrix are carried out in SMIS1. The other operations are 
carried out in subroutines which are called from SMIS1. Most 


of the error messages are called from SMIS1. 
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Figure 1. Functional Flow Diagram 
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Si. Listing of a Matrix (MLIST) 


MLIST stores each matrix by name, number of rows and 
columns and starting position in the array, A. It checks the 
name of a new matrix to be stored against the names of pre- 
viously stored matrices to see if deo mene has been used before. 
If there is another matrix with the same name, the previously 
stored matrix will be deleted before the new matrix is stored. 

4. Printing of Error Messages (ERRER) 

When an error is made in an operational command,.ERRER 
is called. It prints out an error message, deletes the current 
Operational command and returns control to the typewriter. 

ae Locating a Matrix (MFIND) 

Locates the matrix or matrices to be used in an 
operation. It finds the number of rows and columns of a matrix 
and the starting position in the array A. 

6. Calculating Eigenvalues and Eigenvectors (JACVAT) 

Determines all of the eigenvalues and eigenvectors of 
a matrix. It uses the threshold cyclic Jacobi rotation method 
feet. 4). 

He Formation of Matrix N (P) 

The bicubic polynomial, 


p=c x3y34c x3y24c x3y +c x34c x2y34c x2y24c x2y +c x24 
1 2 3 4 5 6 7 8 
c xy3+c. xy24+c xy+c x4c y34c. 0y24+c—Oo¥se 
9 10 11 he ee 14 15 6 


is used in the polynomial fit operation. This polynomial 


can be written as 


ba 


Perey) AC) 
ge a Os eT 


where (F,) = ( x3y3 x3y?_ x3y,x3,x2v3 x2v?_ x2y x2, 
xy3,xy2,xy,x,v3,¥2,y,l ) 


and {C,} = eh lea ane ) aS In the POLFIT operation, values 


C16 
of P at all of the sixteen nodal points of the grid shown in 
Figure 2 must be used to solve for all of the unknown co- 


efficients of the polynomial using the following matrix 


equation: 
eet =en [N} (Cc 
= so eS Lee ae of 
where 
(r(l)) ; 
a: 
(pl2)) 
[IN] = ‘i 
( p(16)) 
al 


Subroutine P forms matrix N for use in the POLFIT operation. 
The bicubic polynomial, P, will be referred to as CUPOL in the 


following sections. 
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Figure 2.° Grid For Subroutine P 
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Bs OPERATION SUBROUTINES 
a Inverting a Symmetric Matrix (SYMINV) 

SYMINV inverts a symmetric matrix, A, by Gauss 
elimination [Ref. 4]. If the matrix is not symmetric, it takes 
the average of the two numbers, A(i,j),A(j,1i), and replaces 
both by the average value and prints a warning that the matrix 
was not symmetric. 

(ae Adding, Removing and Storing a Submatrix (STOSM) 

STOSM is called to do three operations. 

a) Store a submatrix at an indicated position ina 
matrix when called by the command, STOSM. 

b) The command, RMVSM, calls STOSM to remove a sub- 
matrix from the original matrix. 

c) ADDSM calls STOSM to add a submatrix to a matrix 
at an indicated position. 

Sis Function Generation (FUNGN) 
This subroutine takes a function described by a set 


of points (Xi,¥;) , not necessarily equally spaced, and 
PEES2 8 9eon 


calculates, by straight line interpolation, ordinates at equally 
spaced values OL x 
4. Printing of Graphs (PLOT) 
PLOT prints a 10 inch by 10 inch graph of 1 to 4 
curves. It automatically scales the graph to the maximum and 


minimum values of X and Y. The axeS are plotted as 
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with the characters, X and Y, appearing at every inch on their 
respective axes. When all of the points are plotted on one side 
of either axis, an offset axis replaces the real axis and is 
indicated by 

ayes ete A offset X axis 

b) B------ Bo GEfESeEENY axis 
The origin of the axes and the X and Y scale factors are cal- 
culated and printed below the graph. 


Sie Determining the Coefficients of a Bicubic Polynomial 
POLFIT) 


The 16 unknown coefficients of a bicubic polynomial 
surface can be determined in terms of 16 known values of the 
polynomial at 16 known values of the independent variables. 


POLFIT solves the matrix equation, 


iN] 1c} .= [P] 


for the matrix of unknown coefficients, C, where matrix N is 
formed by the subroutine P. Matrix P is a matrix of the values 
Ge CUPOL at each nodal point and is the input to POLFIT. 

6. Printing Contour Maps (CONTUR) 

CONTUR prints a contour map of CUPOL. There are 21 
possible contour lines which are represented by the alphabetic 
feccers, A through J, 0, Q through Z. The input is either the 
coefficients of CUPOL with the upper and lower limits of X and 
meor the values of the function at the grid points used for 
EOLFIT. In the latter case, POLFIT is called and the coeffi- 


Gients of CUPOL are calculated. 
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CUPOL is then scanned over the ranges of X and Y to 
determine the maximum and minimum Z ordinate. The range of Z 
is divided into 42 equal parts represented by alphabetic 
characters and blanks assigned to alternating contour values. 
CUPOL is then evaluated at each point on a line, the correct 
character inserted and the line printed. 


ie Formation of a 2 by 2 Stiffness Matrix (FORMK) 
FORMK forms the stiffness matrix [Ref. 9], 


= 


for a simply supported beam having length, L, Young's Modulus, 


E, moment of inertia, I, shear modulus, G, and effective shear 


area, A. 
The stiffness coefficients are: 
241 2+8 B! 
Ky L ( #5) = k35 
i 2EEe/ t= a: 
=~ BE ( | me Noa 
where 
E 
p = St 
LAG 
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The matrix equation, 


relates end moments, M), Mo, to their corresponding rotations, 


O18 8. as shown in Figure 3. 


Figure 3. Simple Beam 


8. Numerical Integration (ITGRAL) 


The subroutine, ITGRAL, is capable of performing 
numerical integration of one and two dimensional functions 
described by discrete sets of points, either by Simpson's rule, 


or, if the function to be integrated is known, by Gaussian quad- 


rature. 
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oF Numerical Integration of a Set of Differential 
Equations (RESPON) 


In the study of vibrations of systems with many 
degrees of freedom, the following system of ordinary differential 
equations governs the problem: 
[M] {X} + [Cc] {x} + [Kk] {x} = {p} (1) 
where 
[M] is the mass Maer 32, 


[C] is the damping matrix 
and 


[K] is the stiffness matrix. 
In. general, i£ 

{x} = [v] {x} 3 (2) 
where matrix [V] is the matrix of undamped mode shapes and {x} 
is the vector of the time dependent amplitude function for each 
mode, we get after substitution of (2) into Gloue 

mtv) ixhe+ [C) [Vv] {x}. 4+ [RK [vl 4x} = {2} (3) 
After premultiplication by [v]7? and the use of the orthogonality 
properties of normal modes: 


[V]? [mM] [Vv] = [*my 


[V]T [kK] [Vv] = [-k~] = [~w2m_] 


[v]? [c] [v] = [~c_] = [—2dum_] (2) 


(@) such an orthogonality condition imposes some 
restrictions on the form of the damping matrix [C]. 
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where 


2mj 5 


and using the notation 


[v]? {Pp} = {p} 


we get 

[~m.] {x} + [-2Awm_] {x} + [-w2m_] {x} = {p} (4) 
or 

{x} + [-2\0_] {x} + [-w2_]) {x} = {p*} (5) 


where p* = Mas 
PS Ba/ ; 


This last set of equations is completely decoupled 
and can easily be numerically integrated. The subroutine 


RESPON performs this task. 
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III. OPERATIONAL COMMANDS 

The following is a list of the operations that IMIS con- 
tains. The operations must conform to the standard format 
shown in Figure 4. OPERATION stands for any of the symbolic 
command names mentioned below. A, B, C and D are matrices used 
as arguments of the operation. Nl, N2, N3 and N4 are four 
integer arguments that must always be right justified. Sl and 
S2 are two floating point arguments. The arguments can be 
input or output values and are not all necessary for each 
command. 

There are two versions of the matrix interpretive system 
in use at the Naval Postgraduate School. One is for time 
sharing and one for batch processing. The commands denoted by 
w=" Can only be used in conjunction with the time sharing 
option. The READ, WRITE and PUNCH commands make use of 
sequential input/output files which are available with the 
FORTRAN IV language [Ref. 8]. (P) 

The READ command uses the file FTO4F001. This file is 
used to load matrices into IMIS from punched cards and is read 
into the computer prior to executing IMIS. Any data stored in 
the file is automatically erased when new data is read into the 
file. 

File FTO8F001 is used to hae ane output which the user 


wants to print offline. This file is not automatically erased 


(b)A certain knowledge of the CP/CMS time sharing system 
is needed when using a remote terminal of the computer. It is 
presumed that the reader is already familiar with this system 
Peer. 10). 
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when new data is read into the file. New data is appended to 
the data that is already stored in the file. Therefore, the 
file must be erased before executing IMIS. To erase the file, 
the command, 

ERASE, FILE) FTO8F001 

is used. After exiting from IMIS, the command, 


OFFLINE, PRINTCC, FILE) FTO8FO01 


will cause the file to be printed on the offline printer. 

The command, PUNCH, uses file FTO7FO01 to store data to 
be punched onto cards. This file appends new data at the end 
of old data so the file must be erased before using IMIS. The 
command, 


ERASE) FILE, FTO7FOOL 


erases the file. To offline punch the data stored in the file, 
the command, 

OFFLINE,PUNCH)FILE,FTO7F001 

is issued after exiting from IMIS. Examples of the use of 
these files are given in the sample problems. The operational 
commands are: 

START This command causes the elimination of 
all matrices which are in core storage 
and must be the first command used in 
the solution of a problem. if an error 
is found during execution of a command, 


an error message is printed on the 
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typewriter, the current command 
eliminated and control returned to the 
typewriter. If an error is found 
during the execution of a problem with 
DSMIS, the current problem is auto- 
matically deleted and the remaining 
cards scanned until a new START or 
STOP command is found. 

LOAD Matrix Load: 
A = Name of matrix to be loaded. 


Nl 


Il 


Number of rows of matrix A. 


N2 = Number of columns of matrix A. 


N3 Format) controls andicator) fori the 
data. 

The elements of the matrix are 

typed or punched on cards row-wise 

after the LOAD command according 

to one of the options determined 

by the value of N3: 

aera No= I or blank, the format’ vs 
(12F6.0). A maximum of twelve 
numbers per line, using as many 
lines as required to write the 
entire matrix must be used. If 
the decimal point is omitted, it is 
assumed to be at the extreme right 


of the field. 
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PRINT 


STOP 


ZERO 


if N3= 2, (6F12.0), 6 numbers per 
Pine “in -fxelds*o£ “I2 "columns. 
Ti N3S= 3, “(4F 160), 4 numbers per 
line in fields of 24 columns. 
if N3= 4, (3F24.0), 3 numbers per 


line in fields of 24 columns. 


Matrix Print: 


Nl= 


Name of matrix to be printed. 
Number of lines, which follow the 
PRINT command, to be printed as a 
Mabel “for the* matrix. “Column iS 
of the command card is used for 
Carriage control in DSMIS only. 

A one in column 13 will cause the 
matrix to be printed on a separate 


output page. 


This operation causes the number of 


lines designated by Nl, which follow 


the 


operation, to be printed as 


remarks. 


This operation ends a run of problems 


and 


is the last command to be issued 


in a ‘stack Of probilems. “Tt *can be 


used only once per run. 


Formatron “of a Nuri Matrix: 


A 


Name of null matrix. 
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Nl= Number of rows of matrix A. 
N2= Number of columns of matrix A. 

ADD Matrix Addition - [A] + [B] = [A] 
Matrix A is replaced by the sum of 
matrices A and B. 

SUB MatrixaSubiractione— i [A] - [B] = [A] 
Matrix A is replaced matrix A minus 
matrix: B. 

MULT Matrix Multiplication - [A] + {[B] = [Cc] 
The product of matrix A times matrix 
B is generated and defined by matrix C. 

TRANS Matrix Transpose - [B] = [a]? 
The transpose of matrix A is generated 


and defined by matrix B. 


¢ 


INVERT Matrix Inversion - [A] = taj 71 
Matrix A is replaced by the inverse 
Of atsei¢. 

SYMINV Symmetric Matrix Inversion - [A] = zi 
Symmetric matrix A is replaced by the 
inverse of itself. 

STOSM Storey Submatxi x win) a: Matrisc: 

A = Name of large matrix. 

B = Name of submatrix to be stored. 

Nl= Row number in large matrix of 
first element of submatrix. 


N2= Column number in large matrix of 


first element of submatrix. 


Zo 


Elements of A are replaced by 
elements of B in designated area. 
STODG Store Row Matrix on Diagonal of Matrix: 
A = Name of square matrix. 


B Name of row matrix to be stored on 


diagonal of matrix A. 
Diagonal elements of A are replaced by 
elements of B. 
SCALE Scalar Multiplication = [A] = Sil |[A] 
ADDSM Add Submatrix to a Large Matrix: 
A = Name of large matrix. 
B = Name of submatrix to be added 
into large matrix. 
Nl= Row number in large matrix of 
first element of submatrix. 
N2= Column number in large matrix of 
first element of submatrix. 
RMVSM Extracts Submatrix [B] from [A]: 
A = Name of large matrix. 
B = Name of submatrix defined as matrix. 
Nl= Row number in large matrix of 
first element of submatrix. 
N2= Column number in large matrix of 
first element of submatrix. 
N3= Number of rows in submatrix. 
N4= Number of columns in submatrix. 
Matrix A is not modified. 
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LOG 


MSCALE 


RMVDG 


EIGEN 


The Log of each Element 

Each element in the matrix A is 
replaced by, thé natural log of Ehe 
element. 

Scalar Multiplication 

Each element in the matrix A is 
replaced by B times the element where 
B has previously been defined as a 
Al eee al Matrix. 

Extracts Row Matrix B from Diagonal 
©f Square matrix, A. 

A = Name of Square matrix. 


B 


Name of row matrix composed of the 


Giaconal elements of matrix A. 


* Matrix A is not modified. 


Eigenvalues and Eigenvectors: 

The eigenvalues and eigenvectors of 
the system, [A] [X]=A[B] [X], are 
determined where [A] is a symmetrical 
matrix and [B] is a diagonal matrix 
of positive elements stored as a row 
matrix. 

A = Name of matrix A. 


B Name of matrix B. 


C = Name of matrix of eigenvectors 


stored row-wise. 
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SQREL 


INVEL 


DELETE 


Name of row matrix of eigenvalues, 
Xe 

The number of eigenvectors to be 
calculated. The ordering of 
eigenvectors and eigenvalues is 
determined by the sign of Nl as 
follows: 

If Nl is positive, eigenvalues are 
arranged in descending order. 

If Nl is negative, eigenvalues are 
arranged in ascending order. 
Eigenvector normalizing. 

if NZ is) blank: 

Eigenvectors are normalized so 
that [Clip] ic! = (2! 

If NZ is i: 

Eigenvalues are normalized so 


that [A] = [C]~ fa] Ic] 


Square Root of Each Element 


Each element in the matrix A is 


replaced by the square root of itself. 


Inversion of Each Element 


Each element in the matrix A is 


replaced by the reciprocal of itself. 


Matrix Deletion 


The matrix A is eliminated from core 


storage. 
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DUPL 


ENVEL 


RESPON 


Matrix Duplication: 


ha 


B 


Name of matrix to be duplicated. 
Name of matrix defined to be 


identical with matrix A. 


Envelope Value of Matrix 


The 


row 


maximum absolute value in each 

of matrix A is printed along with 
column number times an interval. 
Name of matrix A. 


Interval between columns. 


Numerical Integration of a Set of 


Differential Equations: 


{%}42) [-ur] {x} + bw 2] {x}={s}P (t) 


i = 


with initial conditions, x(0)=x(0)=0. 


Name of a row matrix containing 
the diagonal elements of matrix 
[-—w.].. The dimension of A must 

be (lxn) where n is the number of 
equations in the set. 

Name of a column makeas containing 
the elements of {S}. The dimension 
of B must be. (nx1). 

Name of a row matrix containing 
the values of the function P(t) 
evaluated at equal intervals, At: 


PO), P (AU) P (2K) |. «2 ,P (mat) 
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The dimension of C depends on the 
choice of At which must be small 
for accurate results. The size of 
At should be determined by actual 
numerical experimentation with 

the command. 

D = Name of the output matrix contain- 
ing the calculated values of the 
variables as shown in Figure 5. 

Nl= Output interval i. The values of 
the n unknowns, {x}, are printed 
at 2At, 2ZPAC Sint uc, 2h aa 
all of the values are printed. 

Sl= Value of the damping ratio, i. 

~$2= Time increment, At, used in the 


numerical integration. 


x4 (1At) ,x, (21At) ,x, (31At), oienoeiie 


Xo (iAt) , x5 (2iAt) ,x, (3idt), aces 


e 
e 


¢ 
e 


x, (1At) ,x, (21At) ,x, (3iAt), 50 34 


Fig. 5. OUTPUT MATRIX, D, FOR RESPON 
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FUNGN 


FORMK 


Function Generation: 
A = Name of matrix which defines a 


function in terms of line segments 


B = Name of matrix to be formed which 
is composed of the y ordinates at 
equal x ordinate intervals of the 
function defined by matrix A. 

Nl= Number of y ordinates to be 
generated. The first ordinate 
will be equal to Ya: 

Sl= X ordinate interval. 


Forms a 2x2 Element Stiffness Matrix: 


A = Name of matrix to be formed. 


The line following the command is as 

follows: 

a. Columns 1-12 Moment of inertia of 
member, I. 

b. Columns 13-24 Effective shear 
area of member, A. 

e. Columns 25-36 Length tof member, Ll. 

d. Columns 37-48 Modulus of 
elasticity of member, E. 

The shear modulus, G, is assumed to be 


equal to E/2.4. 


Sill 


* 


* 


* 


PUNCH 


READ 


CONT 


WRITE 


Punches Matrix onto Cards: 

A = Name of matrix to be punched. 

Nl= Number of comment cards to be 
punched as a heading for matrix A. 

After the PUNCH command is issued, the 

next line contains the format elected 

by the user as follows: 

(NNTRR.SS) 

where 

NN= Number of elements per card. 

RR= Field width per element. 

SS= Number of significant digits. 

T = Type of data field. T must be a 
standard F or D Fortran format. 

Loads Matrices from file FTO4F001 

The load operation card must have a 

I. in Si. sthe dast card musi bela 

Continue Operation card. 


Continue Operation 


Indicates that all matrices in file 
FTO4F001 have been loaded into IMIS. 
Writes Matrix into File FTO8FO001: 
A = Name of matrix to be written. 


Nl= Number of lines of remarks to be 
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SoS minty — A 


PLOT 


N2 

Plot 
Sets 
A, B 


a 


Matr 
have 
have 
with 
ates 
ord1 
one 


may 


written as a label for the matrix 

A. 

must be equal to l. 

s One to Four Curves Described by 
of Points Contained in Matrices 
q2tGiyaeDs 

First matrix to be plotted on one 
graph. ‘Plotted asia, "+". 

Second matrix to be plotted on 

one graph.  Pillotted asa. ™*™. 

Third matrix to be plotted on one 
graph... Plotted as a "O".; 

Fourth matrix to be plotted on one 

graph. Plotted as a".". 

lL. “Offline, prink indicator. Graph 

will be stored in file FTO8F001 

Number of lines of remarks to be 

printed at bottom of graph. 

ices to be plotted do not have to 
the same number of points. They 
to be entered as N by 2 matrices 
the first column the X co-ordin- 
and the second column the Y co- 
nates. When plotting more than 
curve on the same graph, a conflict 


arise if more than one curve passes 
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ie ea 


through the same point. The symbol 

printed at that point is determined by 

the following. order: D,.C, B and A; 
POLFIT Determines the Coefficients of CUPOL: 


A = input, matrix. (sx) 


B = Output. matrix,of, coefficients of 
CUPOL, (Cy,Co,--+sCy¢) 
where 
A(1) = H, 
2) =H 
A(2) y 
Le — A4 —— AS ——A6—.—} X 
Hy a 
n aa eins eae 
All Pee a satis 
| Tr ALO eH) ae 
NE 
and 


3-3 au 3 3 
CUPOL=C)X Y +CyX Y +C3X Y+C4xX + 


3 22 2 2 
CuX7¥ 340, X2¥ 240, x2v4+C 4x24 
C.xXY¥34+C, .xy24+C. .XY+C. X+ 

9 10 1l 12 


¥3+¢C_, ,Y24+¢C 


“13 14 15 


Y+Ci¢ 
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CONTUR Contour “Maps “CUPOL: 
N2= Number of lines of comments. 
N3= Input type indicator. 
a. NG= Blank or 1. 
A= TApUE matrix. Same as for 
POLFIT command. (18x11) 
B = Output matrix of coefficients 
of CUPOLE. “(16xr) 
b. N3= 2 
A= Input’ matrix. (20x1) 
A(1) = Lower bound of X. 


A(2) 


Upper bound of X. 


A(3) Lower bound of Y. 
A(4) = Upper bound of Y. 
A(5) through A(20) are the 
coefficients of CUPOL. 
* N4= Offline print indicator. Stores 
contour map in file FTO9FOO1. 
ITGRAL Integration of Functions Represented 
by Discrete Values or Polynomials: 
A = Input matrix. (Rx1) 
N3= Integration type indicator: 
a. If N3= 1, Simpson's Rule for one 
dimension. R must be even. 
A(i) = H, 
A(2) through A(R) are evenly 


spaced values of function. 
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If N3 = 2, Simpson's Rule for two 


dimensions. 
A(1) = A, 
A) = Hy 


A(3) through A(R) are values of 
function entered row-wise as for 
POLPEET. 
R = N2+2, N must be odd. 
There must be the same number of 


values in both dixections. H. and 


x 
Hy can be unequal. 

If N3= 3, Gauss Quadrature for one 
dimension. 

A(1) = Lower bound of X. 

A(2) = Upper bound of X. 

A(3) through A(14) are the co- 


efficients of an llth degree 


polynomial: 


Ee 10 9 8 d/ 6 
AX tA,X tA X tAEX tADX tAgX + 


2 
D Saale Bi. X+Aj 4 


5 4 3 
AgX tA, 9% +A, 4% tA, 9 Tye 


If N3 = 4, Gauss Quadrature for two 


dimensions. 
A(1) = Lower bound of X. 
A(2) = Upper bound of X. 


A(3) = Lower bound of Y. 
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Ce 


‘7 to 


Bkorei 


A(4) = Upper bound of Y. 
A(5) through A(20) are coefficients 


of CUPOL. 
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IV. LIMITATIONS 


A. SIZE AND NUMBER OF MATRICES 

The storage available for all of the matrices in any 
problem is 80,000 bytes or 10,000 double words. This limit 
is arbitrary and controlled only by the total core storage 
available at any installation. For the naval Postgraduate 
School installation, 10,000 double words was selected in order 
to reduce turn around time in batch processing and reduce the 
size of the data set required for time sharing. The total 


number of matrices which can be stored is 100. 


B. MATRIX NAMES 

Each matrix is referred to by its name composed of one to 
six alphameric characters including blanks. The name is select- 
ed by the user when the matrix is first loaded or generated by 
an operation. If the name for a previously defined matrix is 
used later as a name of a new matrix, the former matrix is 
deleted before the operation is executed and is replaced by the 
new matrix. This is done so that there is only one matrix with 


a certain name at any one time. 


es PLOT 

There are 200 points in the Y direction on the graph. The 
X direction has 130 points for the remote terminal printing and 
160 points for the offline printer. The difference is due to 
the different number of lines per inch. The accuracy of the 
position of a point depends on the range in both the X and Y 
G@irection. The maximum error in the Y direction is the range 


of Y divided by 200. The maximum error in the X direction is 
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the range of X divided by either 130 for the remote terminal or 
M60 for the offline printer. If the X and Y ordinates of two 
or more points fall on the same nodal point on the graph, only 
one point will be plotted. If the points are from different 
Matrices, the matrix hierarchy (see page 34) determines which 


will be plotted. 


Dd. INTEGRATION 

Integration by Simpson's Rule is restricted to equally 
Spaced points in any one direction. The error is of the order 
of H2 for integration in one dimension and of the order of H2 
times HE for two dimensions. H, does not have to equal Hy but 
there has to be an odd number of points in any one direction. 
Integration by Gauss Quadrature for six Gauss ordinates is 


exact for all polynomials of degree eleven or less in one dimen- 


sion and for CUPOL in two dimensions [Ref. 11] 


E. POLFIT AND CONTUR 

The polynomial, CUPOL, will pass through the values of 
the function at the nodal points of the grid but may not be a 
true representation in the space between nodal points. If the 
range of values of the function is less than 10-2, a contour 


map will not be printed. 
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V. HOW TO CALL THE MATRIX INTERPRETIVE SYSTEM 


A. TIME SHARING (IMIS VERSION) 

At the Naval Postgraduate School, IMIS is stored ina 
library (AED) with the name, LIT, so that any number of users 
can use the system at any one time. In order to call LIT, it 
must be called as a subroutine. A Fortran program called IMIS 
has to be written as follows: 

CALL LIT 

STOP 

END 

There are two other libraries which must be loaded along 
with AED library in order to use LIT. This can be done with an 
executive program called SMIS as follows: 

FORTRAN IMIS 

GLOBAL LOADER TXTLIB AEDLIB SYSLIB 

LOAD IMIS 

START 

SMIS is then executed with the command, $ SMIS. The 
computer will then type out the executive program with each 
command preceded by a time. Next it will type 

EXECUTION BEGINS ... 

Ac this time you are in IMIS and ready to begin the first 


problem. The first command will be START. 


B. BATCH PROCESSING (DSMIS VERSION) 
The first card must be a job card for the batch processing 


system as follows: 
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Wace Sela: a! 


ip POS 6 ort oe i 


//aaaznnmp,JOB(nnnn,SSSssFP,TTTT) , "IDENTIFICATION' 
where 


aaa = The £irst three. letters of the user's Last name. 


Z Job sequence number 

nnnn= User's number. 

Ssss= Project number. 

TTTT= Student section or faculty number. 

The next six cards are control cards. Card five is used 
only when using the punch command. The cards are as follows 
stating in column one: 


//JOBLIB,DD DSNAME=DSMIS.F59CI,UNIT=2314,DISP=(OLD,KEEP) , 


b 


//ppbbbbbbbbbbbYOLUME=SER=MARY 
//ppEXEC)PGM=DSMIS , REGION=160K 
//FT06F001,DD,SYSOUT=A 
//FT07F001;,DD,SYSOUT=B 
//FTO5F001pDDp* 


The ag is where blanks are left. 


The problem deck follows with the first command a START 
command and the last a STOP command. The last card of the 
deck is a standard delimiter card: 

/* 

On the job request card, check the item marked DISK PACK 

under the SPECIAL INPUT heading and write in "MARY". If the 


punch command is being used, under the heading SPECIAL OUTPUT 


check the item marked PUNCHED CARDS, 
Ani v 


( Loi Sis) 
| 


VI. SAMPLE PROBLEMS 


A. THE EIGEN SYSTEM 


The matrix equation 


mm) {XxX} = AB 1x} 


where 
0:0 -1.0 0.0 
fA} = -1.0 20 =1 0 
0.0 -1.0 20 
and 
i. 5 0.0 0.0 
fei = 050 2.0 0.0 
0.0 OR0 225 


is to be solved. The entire sequence of events is as shown 
below. Everything written in lower case was typed by the 


operator. The computer response is in upper case. 
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B. CONTOUR MAP 
A contour map for a surface described by the ordinates at 


the sixteen nodal points of the grid shown on page 47 is 


desired. 
10'.0 1250 8.0 LOY. 
PA} = 2.0 LS 97.10 8.0 
8.0 9.0 Pe. 0 1230 
£O).:0 8.0 P2210 LOGO 
and 
H,=Hy=1 


The results are as follows: 
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The command 
offline punch file f£t07f001 
will cause the matrix B to be punched onto cards with the 
indicated format of four numbers per card with a field width 
@= 15 and 8 significant figures. 
If the input for the contour map had been the coefficients 


of CUPOL, the first page of output would have been as follows: 


THIS IS THE FIRST PAGE OUTPUT COEFFICIENTS AS INPUT. 


=0.333300 X3Y3 1.500000 X3Y2 =2.833000 X3Y 2.000000 X3 
1.500000 X2Y3 -6.750000 X2Y2 12.750000 X2Y -9.000000 x2 
=2.833000 XY3 12.750000 XY2 -18.750000 XY 9.000000 xX 
2.009000 Y3 -9.000000 Y2 9.000000 Y 10.000000 
Q= 7.083 7699378 RANGE( 0.708377D Ol 0.1295°94D 02 ) 
A= 7.3776 B= 7.6713 C= 7.965) O= 8.2589 E= 8.£527 
F= 8.8465 G= 9.1402 H= 9.4340 T= 9.7278 J= 10.0216 
Q= 10.3154 R= 10.60°2 $= 10.9029 T= 11.1967 U= 11.4905 
v= 11.7843 W= 12.0781 X= 12.3719 Y= 12.6656 z= 12.9594 
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Cc. GRAPH 


BM plot of the: chunction 


x = e St sin(54.99t) 
is desired from t = 0.0 to t = 0.3. The values were calculated 
in a separate program and the values of X were punched onto 
cards. They were prepared for offline reading into file 
FTO4F001. 
The first card is a standard CP67 offline read card as 
follows starting in column one. 


CP67USERID)) GUxxSSSS))))ppTTTT 


where 


x Terminal number user will be using. 


ssss= User's account number. 
TTTT= User's name. 
The second card is as follows: 


OFFLINE}, READ) FILE,FT04F001 


the next card is a LOAD command card with a 1. in Sil. 
LOAD A 80 2 5 

The matrix is punched onto cards row-wise with the same 
format as one of the options under IMIS and these cards follow 
the LOAD card. As many matrices as desired may be loaded in 
this manner. The last card is a continue card as follows: 
CONT 1 

This card indicates that all matrices have been read 


from file FTO4F001. The deck is then loaded into the computer. 


Si 


After logging onto the time sharing system, the command to read 
the file is as follows: 
offline read * 

OFFLINE READ READ FILE FTO4F001 


IMIS is then executed and the sequence is as follows: 
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The command 
OFFLINE PRINTCC FILE FTO8FOOL 


will cause the matrix A to be printed on the offline printer. 
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VIL. | IMPROVEMENTS 


This system is an open ended system, i.e., any matrix 


operation with fewer than four matrix arguments, four integer 


arguments and two 
command and added 
sequential. Each 


no facilities for 


scalar arguments can be coded as a new 
to the system. As coded this program is 
command is executed once after it is read, 


automatic repetitive execution of one or more 


commands is included. Such a feature would be very useful if 


it could be implemented without introducing extraneous 


complexities to the language. 
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